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ABSTRACT 

In this Letter, we present a model connecting the cosmic star formation rate (CSFR) to 
the growth of supermassive black holes. Considering that the evolution of the massive 
black hole is dominated by accretion (Soltan's argument) and that the accretion 
process can be described by a probabilistic function directly regulated by the CSFR, we 
obtain the evolution of the black hole mass density. Then using the quasar luminosity 
function, we determine both the functional form of the radiative efficiency and the 
evolution of the quasar duty-cycle as functions of the redshift. We analyze four different 
CSFRs showing that the quasar duty-cycle, S(z), peaks at z ~ 8.5 — 11 and so 
within the window associated with the reionization of the Universe. In particular, 
^max ~ 0.09 — 0.22 depending on the CSFR. The mean radiative efficiency, fj(z), peaks 
at z ~ 0.1 — 1.3 with 77 max ~ 0.10 — 0.46 depending on the specific CSFR used. Our 
results also show that is not necessary a supercritical Eddington accretion regime to 
produce the growth of the black hole seeds. The present scenario is consistent with 
the formation of black hole seeds ~ 1O 3 M at z ~ 20. 
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nuclei — quasars: general 



1 INTRODUCTION 

There is strong evidence that nearly all galaxies contain 
supermassive black holes (SMBHs) in their centers, and that 
the evolution of the supermassive black hole and its host 



galaxy are connected (see, e.g., Ferrarese & Merrit 2000). 
In particular, the masses of the black holes range from 
~ 10 6 Mq for galaxies with small bulges up to ~ 1O 9 M0 for 
galaxies in cores of groups and clusters of galaxies (see, e.g., 
Margorrian et al.|1998[ ). These SMBHs are present not only 
in the local Universe, in the form of the low luminosity active 
galactic nuclei (AGN), but also in the early stages of galaxy 
formation as can be seen from quasars discovered beyond 



z > 6 (see, e.g., Fan et al. 20031. Furthermore, accretion 
onto massive black holes is generally accepted as a way to 
power strong emission as observed in the cases of AGNs and 
quasars. 

On the other hand, it is reasonable to consider that the 
growth of SMBHs can be regulated, in some way, by the 



cosmic star formation rate - CSFR (see, e.g., 


Franceschini 


etaLj [19991 |Haiman et al.||2004| [Heckman et al. 2004; 


Merloni et al. 


|2004| |Mahmood et al. 


|2005| |Wang et al.| 



|2006[ ). This is reinforced by the fact that the CSFR can be 
directly connected to the masses of dark matter halos (see, 
e.g., IConroy & Wechsler||2009| iPereira & Mirandal 120101). 



The dark halos are the natural nursery for the birth and 
growth of the massive black holes. In principle, measuring 
the masses and accretion rate of the black holes that 
drive the AGNs and quasars could help us understand the 
evolution of these sources, their connection with the CSFR, 
and the contribution of mini-quasars to the reionization of 
the universe. Here, we present a formalism permitting to 
confront several CSFRs with the quasar luminosity function 
(QLF). This formalism could also be used for a better 
estimate of the CSFR up to redshift ~ 7 using the QLF 
as an observational data to be fitted by the 'theoretical' 
CSFRs. Furthermore, this work could contribute to the 
study of the feedback processes associated with both: star 
formation at higher redshifts and growth of supermassive 
black holes. As main results, we derive the functional form 
of the radiative efficiency associate to the accretion process 
of these black holes and the evolution of the quasar duty- 
cycle. Through this paper we consider standard cosmological 
model (ACDM) with f2 b = 0.04, fi m = 0.27, fi A = 0.73, 
h = 0.73. 
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2 SMBHS AND THE CSFR 



We consider that black holes grow by accreting matter 
(Soltan's argument; see, e.g., Soltan 1982 Wang et al. 
|2008[ ) and that the accretion process can be described by a 
probabilistic function directly regulated by the CSFR. This 
can be described in the following way 



/ \ o P*bh(z) 

PBh(z) = PbH" 7- TyT, 

P*bh(2 = 0) 



(1) 



where p BH represents the black hole mass density in our 
local universe. The function p*bh(z) makes the connection 
between the CSFR and the black hole mass density at 
redshift z. In particular, we consider 



P*bh(z) = 



dz' 



(2) 



In Eq. pj, p*(z) represents the CSFR at redshift z, 
P{td) is the probability per unit of time of the black hole 
grow up by accreting matter from the environment, and the 
(1 + z) term in the denominator considers the time dilatation 
due to the cosmic expansion. The time delay td makes the 
connection between the redshift zt at which the accretion 
disk forms around the black hole and the redshift z at which 
the material is incorporated into the black hole. Thus, we 
have 



td 



9.78/i _1 Gyr 



(i + z')V n A + n m (i + «') 3 



dz' . 



(3) 



We consider that P(td) = Bt^ and this function 
is normalized as J t * now P(i d )cfe d = 1, where i m i n is the 
minimum time required for the accretion process to take 
place and t nov is the age of the universe. We assume 
tmin = t s = 4.2 x 10 7 yr, where t B is the Salpeter time. 
It is worth stressing that P(td) oc t d has been used by 
different authors in several astrophysical contexts (see, e.g., 



Rcgi mbau fc de Freitas P achcco 2006; Rcgi mbau fe Hughes| 



2009 and references therein). 



In this work, we consider the CSFR, p*(z), derived by: 
PM - |Pereira fc Miranda] ( |2010[ > who obtained this function 
from the hierarchical scenario for structure formation 
using a Press-Schechter-like formalism; SH - |Springel fc| 
|Hernquist| ( 2003 ) who obtain the CSFR from hydrodynamic 



simulations; HB - Hopkins fc Beacom| (20061 and F 



Fardal et al. ( 2007 I who derived the CSFR using the 



available observational data. In Fig. [T] we summarize all 
these different CSFRs. 



In Fig. [2] can be seen the black hole mass density as 
a function of the redshift. The results are presented for 
four different CSFRs, as explained above, and using the 
formulation contained in Eqs. |TJ-([2|. Observe that for the 
CSFR derived by PM, pbh/Pbh peaks at z = 1.1 (1.5) for 
n = —1.0 (n = —1.5). Considering n = —1.0, note that PM- 
CSFR produces higher values for pbh(z) at redshifts z > 4 
than those produced by others CSFRs. We also plot in Fig. 
[2] psH (z) derived from the standard model (SM). In this 
case, Pbh(z) is obtained from the integration of the QLF, 
$b(ib, z), in the following way (see, e.g., Small & Blandford 




Figure 1. Cosmic Star Formation Rate (CSFR) derived by: PM 
- |Pereira fc Miranda| (|2010); SH - |Springel fc Hernquist] (2003) ; 
HB - |Hopkins fc Bea,com| 1 2006| ; F - |Fardal et al.| (2007 1. The 
observational points were taken from |Hopkins | (2004|[2007) . 




Figure 2. The co-moving massive black hole mass density. We 
present pbh( z )/Pbh:( 2 ) usm S f° ur CSFRs presented above. 
The results are shown for two values of the n-exponent of the 
probability function (n = —1.0 and n = —1.5). We also present 
for comparison Pbh derived from the standard model (SM). 



T992l|Yu fc Tremaine 2002; M erloni et al.|2004||Wang, Chen| 
fc Zhang|2006| |Shankar||2009| |Shankar et al.|2009) 



Pbh (2) 



dt 

dz 



Li 



c 2 rj 



L b $ h (L h ,z)dL h , (4) 



where Lb is the bolometric luminosity taken from |Hopkins| 
et al. (20071, rj is the radiative efficiency and c is the speed 
of light. Note that in order to solve equation Q some 
consideration about rj value should be made. In general, 77 
is assumed constant with typical value ~ 0.1. However, we 
will show in the next Section that rj must be a function of 
time if the growth of the SMBHs is regulated by the CSFR. 
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3 THE MEAN RADIATIVE EFFICIENCY AND 
THE QUASAR DUTY-CYCLE 

The bolometric luminosity of a black hole with accretion 
rate M a is given by L = fjM^c 2 (with fj the mean radiative 
efficiency). Considering that M a is related to the mass 
variation of the black hole by M a = Mbh/(1-0 — 77), and 
that fj is a function only of the redshift 2, we obtain: 



U = c z 



n 



n) 



PBH, 



(5) 



with U being the luminosity density. Deriving the equation 
|l]) in 2 and using the result in |5| produces: 



(6) 



The temporal dependence of fj reflects the fact that 
the active supermassive black holes have a finite lifetime 
(Davis & Laor 2011). In order to derive fj(z), we shall use 
the luminosity density U(z) obtained from the integration 
of the QLF in the following way 

U(z) = f L b $ h {L h ,z)dL b , (7) 
where L b is the lower limit of the QLF. Here, we take L£ 



and <!>b(Lb, 2) from Hopkins et al. (2007 1 . Then, defining 



/(*) 



1-fj' 



we can write 



U = c 2 f{z) Pi 



(8) 



(9) 



Now, we define a functional f'(z, n*i) which will be used 
to map f(z) given by Eq. |8|. In particular, n\ is a vector of 
parameters and so 



f'(z,b 1 ,b 2 ,t q ) = Co 



tu{z) 



+ 



tu(z) 



b 2 



(10) 



where Co is a normalization constant which gives fj(z = 0) = 
ff (with ff =0.1, see, e.g., Hopkins et al.|p007| , i q can be 
understood as a characteristic time-scale, h (i = 1,2) are 
dimensionless constants, and t u (z) is the age of the Universe 
at redshift 2. 



The parametric form of Eq. ( 10 1 is widely used in the 
literature. For example, Eq. (12) of Hopkins & Hernquist 
(2009) and Eq. (13) of |Hopkins et al.| ( j2006|) are similar to 
Eq. ( 10 1 presented here. Furthermore, Wang et al. ( 2006 ) 



used the same parametric form to determine the mass 
function of supermassive black holes. Note that, we are using 
the Eq. |7| as a source of information to obtain fj(z). Thus, 
using Eq. | |10[ ) in Eq. |9| is possible to write 



ei(ni,(7i,2i) = \\U(zi) - c f {z h m)pnB.{zi)\\. 
The function that will be minimized is 



(11) 



Table 1. Best fit parameters of r](z). 



n 


CSFR 


61 


b 2 


t q (Gyr) 


-1.5 


PM 


2.85 


2.39 


5.74 


-1.5 


HB 


1.51 


1.29 


13.19 


-1.5 


F 


0.82 


1.47 


5.99 


-1.5 


SH 


2.74 


3.17 


4.93 


-1.0 


PM 


1.81 


1.96 


5.40 


-1.0 


HB 


1.91 


0.74 


14.81 


-1.0 


F 


0.46 


1.06 


4.23 


-1.0 


SH 


2.37 


2.72 


4.72 



Table 2. Maximum values for the mean radiative efficiency. 



n 


CSFR 


'/max 


^max 


-1.5 


PM 


0.38 


1.1 


-1.5 


HB 


0.10 


0.1 


-1.5 


F 


0.11 


0.6 


-1.5 


SH 


0.46 


1.2 


-1.0 


PM 


0.23 


1.1 


-1.0 


HB 


0.10 


0.3 


-1.0 


F 


0.11 


0.7 


-1.0 


SH 


0.39 


1.3 



(12) 



See that fj = f(z)/(1.0 + f(z)) and f(z) = f(z, n ihcst ), 
where nlbest gives the best fit from the equation (12). In 



Table [T] we present the best fit parameters which permit to 
derive the function f](z). 

Figure [3] presents the luminosity density of this work 
when compared to that obtained f rom Eq. ^f\, wh ich comes 
from the integration of the QLF ( jHopkins et al.|[2007[ ). In 
particular, for n = —1.0 we see that PM and SH CSFRs 
produce an excellent agreement with U (2) up to redshift 
~ 6. In the case n = —1.5 we verify that F-CSFR has the 
best agreement with U(z) at 2 ^ 1. On the other hand, from 
2 ~ 1 up to redshift ~ 6, PM and SH, beyond F-CSFR, have 
good agreement with the integrated QLF if we consider the 
parametric form given by Eq. (10). 



In Fig. [4] we present the mean radiative efficiency 77(2) 
as a function of the redshift while in Table [2] we show the 
redshift where fj(z) peaks. See that 77 max is within the range 
0.10 — 0.46 depending on the specific CSFR. In particular, in 
these cases we have 2 max in the range 0.1 — 1.3. These results 
are concordant with accreting black holes which could reach 
fj ~ 0.2 for the most massive systems ( |Narayan|2005[ ) . 

On the other hand, the Eddington mass can be written 
as M c dd = AfeH/is- The bolometric luminosity- weighted by 
the Eddington mass is L = (r/c 2 mMBn) /t B , where m = 
M a /M e dd is the dimensionless accretion rate. Assuming that 
< m > is a function only of the redshift (see |Hopkins et ah] 
2007 1 we can write 



U{z) = 



77 < rh > pbhc 



(13) 
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Table 3. Maximum values for the quasar duty-cycle. 



10 s 



10° 

U 

© 

Sio 7 



10 s 



10 c 



© 

Sio 7 



CSFR 



-1.5 


PM 


0.13 


11.1 


-1.5 


HB 


0.15 


8.5 


-1.5 


F 


0.09 


9.9 


-1.5 


SH 


0.16 


9.4 


-1.0 


PM 


0.19 


10.9 


-1.0 


HB 


0.20 


9.0 


-1.0 


F 


0.14 


9.7 


-1.0 


SH 


0.22 


9.6 



: 1 

: »=-1.0 


111: 




— HP 






— PM 






F 






HB 






SH 






i 



; 1 
: n=-l.B 


iiz 




— HP 






— PM 






F 






HB 






SH 






i 



Figure 3. The quasar luminosity density as a function of the 
redshift. HP represents U(z) obtained from the integration of the 
QLF. Different U(z) directly derived from the CSFRs are also 
presented for n = —1.0 (upper panel) and n = —1.5 (lower panel). 



Using equations ( |13| l and (jSjl, we obtain 
is Pbh 



< to > = 



(1 -77) pbh' 



(14) 



In Figure [5] we present the evolution of < to > with 
the redshift. We can see that the accretion processes are 
more active at higher redshifts for all CSFRs studied in 
the present paper. However, our results also show that the 
accretion processes never reach a super-Eddington regime. 
The accretion process associated to the growth of SMBHs is 




Figure 4. Mean radiative efficiency as a function of the redshift. 
The upper panel presents the results for n = —1.0 while in the 
lower panel we have n = —1.5 for the slope of the probability 
function. 



'well-behaved' during all time associated from the seed black 
hole formation up to the present time. Another point can 



be derived from the equation (14 1. It is possible to find the 



mean accretion- weighted by the lifetime of a SMBH (Merloni 
|et al.|20 04). In particular, we have 



dt , , 
< m > — — az . 
dz' 



(15) 



On the other hand, the quasar duty-cycle, S(z), 
associated to the SMBHs can be derived from the ratio of 
tdc(z) to the Hubble time. This result is presented in Figure 
[6] See that 5(z) is also defined as the fraction of active black 
holes to their total number. As pointed by Wan g, Chen fc| 
|Zhang (20061, this parameter is a key to understand how 
many times and how many black holes are triggered during 
their lifetimes. In Table [3] we show the redshifts where 6(z) 
peaks for each CSFR studied here. Note that typically, S 
is maximum in the redshift range 8.5 — 11.0 and this is 
very curious because the redshift of reionization is 10.5 ± 1.2 
( Jarosik et al.|2011 >. 



4 DISCUSSION 

In this Letter we present a model to obtain the mass 
density of SMBHs from the cosmic star formation rate. 
The key point to do that is to consider Soltan's argument 
and that the accretion process can be described by a 
probabilistic function directly regulated by the CSFR. Our 
model permits to determine the function pbh(z), the mean 
radiative efficiency associated with the growth of the black 
holes, fj(z), and the quasar duty-cycle S(z). In the literature, 
the common way to obtain the mass density of SMBHs 
is by integration of the QLF as presented in Eq. Q and 
considering 77 as a constant. However, here we present a 
different scenario. In particular, we derive pbh(z) from the 
CSFR and then we use the quasar luminosity density in 



Supermassive Black Holes 5 





n = -1.5 














— PM 






HB 






F 






SH 





Figure 5. The mean dimensionless accretion rate for n 
and n = — 1.5 as functions of the rcdshift. 



-1.0 Figure 6. The quasar duty-cycle derived for n 

n = -1.5. 



-1.0 and 



order to obtain the mean radiative efficiency as a function 
of the redshift. From fj(z) is straightforward to obtain the 
dimensionless mass accretion rate and the quasar d uty-cycle. 

If we consider p% K = (5.9ft 3 ) x 10 5 M©Mpc~ 3 ( |Graham 
fc Driver|[2007| |Vika et al.1|2009[), then our model returns, 
using the CSFR of jPereira fc Miranda| ( |2010| ), pgg d = 4.60 x 
10 3 M Q Mpc~ 3 (with n = -1.0). This result is compatible 
with black hole seeds ~ 10 3 Mq at z ~ 20. We also note 
that the quasar duty-cycle has a maximum value close to 
z ~ 8.5 — 11 and so within the observational uncertainties 
associated to the redshift of reionization. As the main 
component of our model is the CSFR, this scenario offers 
several future possibility of investigations. One of them 
is related to the form of the probabilistic function which 
permits to determine the growth of SMBHs from the CSFR. 
Insofar as both star formation and growth of SMBHs can 
be described in the same way, the predictions of different 
probabilistic functions, P(td), can be confronted to different 
observables at high redshifts. Finally, our results also show 
that is not necessary a supercritical Eddington accretion 
regime to produce the growth of the black hole seeds. 
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